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Abstract — We  consider  waveform  design  for  multiple-input, 
multiple-output  radar  systems  for  the  case  where  the  signal,  dur¬ 
ing  propagation,  undergoes  phase  perturbations.  We  formulate 
an  iterative  algorithm  to  obtain  both  waveform  parameters  and 
the  weights  of  the  adaptive  matched  filter.  An  example  of  a  clutter 
and  target  model  is  provided  to  show  how  the  optimal  waveform 
design  improves  the  detection  performance  of  a  random-phase 
radar  compared  to  traditional  waveforms. 

I.  Introduction 

There  has  been  much  recent  interest  in  waveform  design  for 
multiple-input,  multiple-output  (MIMO)  radars.  For  example, 
the  works  in  [1]— [3]  find  the  optimal  waveforms  that  maximize 
the  mutual  information  between  the  received  signal  and  target 
amplitude.  When  the  resolution  capability  of  the  MIMO  radar 
system  is  of  interest,  the  transmit  waveform  can  be  designed 
to  sharpen  the  radar  ambiguity  function  [4],  [5].  However,  the 
obtained  waveform  does  not  necessarily  improve  the  accuracy 
of  target  estimates;  this,  in  turn,  has  motivated  waveform 
design  to  minimize  the  Cramer-Rao  bound  [6],  [7]. 

In  the  context  of  target  detection,  the  signal-to-interference- 
plus-noise  ratio  (SINR)  is  a  well  accepted  design  metric  [8]— 
[11]  for  waveform  design.  In  [8]  the  transmit  waveform  code 
is  assumed  to  be  represented  by  an  unknown  code  matrix 
and  the  code  design  is  formulated  as  an  optimization  problem 
with  several  non-convex  constraints.  A  robust  technique,  using 
minimax  optimization,  was  formulated  in  [9]  for  the  case  of 
unknown  target  parameters,  such  as  Doppler  and  time  delay. 
In  [10],  Friedlander  optimizes  the  values  of  the  discrete  Fourier 
transform  of  the  transmitted  signal  to  maximize  SINR.  Another 
optimization  framework  based  on  SINR  was  proposed  in  [11] 
where  an  iterative  algorithm  is  used  to  estimate  the  optimal 
waveform  as  well  as  the  optimal  receiver  transfer  function.  A 
common  theme  to  all  these  SINR-maximization  approaches  is 
the  need  for  second  order  statistics. 

To  establish  realistic  second-order  statistics,  we  consider  the 
example  of  a  radar  operating  in  a  plasma  medium.  The  medium 
is  assumed  to  contain  plasma  density  inhomogeneities  that 
impart  a  randomization  (denoted  as  scintillation)  of  the  phase 
in  the  spatial  and  temporal  domains.  The  phase  scintillation 
manifests  itself  as  an  apparent  angular  and  Doppler  spreading 
of  the  signal  that  can  be  characterized  by  a  spatial-temporal 
correlation  function.  This  correlation  function  determines  the 
nature  of  the  SINR  and  the  resulting  waveform  and  adaptive 
filter  optimization.  A  key  feature  of  the  plasma  example  is 
that  well-established  fluid  and  electromagnetic  properties  of 


the  plasma  can  be  used  to  determine  the  nature  of  the  spatial- 
temporal  correlation  function  from  first  principles  without  the 
need  to  adopt  empirical  models.  Specifically,  we  consider  a 
class  of  plasma  with  high  Reynolds  number,  which  tends  to 
develop  turbulence  structure  with  a  characteristic  Kolmogorov 
inverse-power  law  wavenumber  spectrum. 

This  paper  proposes  a  joint  waveform  optimization  and 
adaptive  processing  for  a  random-phase  radar  system.  The  sig¬ 
nal  model  for  our  random-phase  radar  system  is  first  discussed 
and  the  SINR  for  the  given  structure  is  derived.  Afterwards,  an 
iterative  optimization  algorithm  is  formulated  to  obtain  both 
the  waveform  and  adaptive  matched  filter  (AMF)  weights.  The 
superiority  of  the  new  joint  framework  to  traditional  waveform 
design  is  then  justified  by  evaluating  the  detection  performance 
of  the  radar  structure. 

The  rest  of  this  summary  is  organized  as  follows.  Section  II 
reviews  the  signal  and  clutter  model.  Based  on  these  models, 
the  SINR  is  derived  in  Section  III  and  the  receive  adaptive 
weights  and  waveform  are  optimized.  Section  IV  then  presents 
numerical  results  that  illustrate  the  gains  possible  using  wave¬ 
form  optimization.  Finally,  Section  V  concludes  the  paper. 

II.  Signal  Model 

The  radar  system  considered  in  this  paper  is  configured 
in  monostatic  mode  with  Nt  transmitters  and  iVK  receivers 
that  are  closely-spaced  in  a  two-dimensional  space.  For  ease 
of  exposition,  the  transmitter  and  receiver  are  assumed  to  be 
configured  as  uniform  linear  arrays  aligned  with  the  x-axis  of 
a  Cartesian  coordinate  system,  with  inter-element  spacings  of 
d,T  and  d r,  respectively. 

The  transmitted  waveform  by  the  m-th  transmitter  is  given 
by 

I 

Um(f)  =  ^2  CrniS{t  -  iTr),  (1) 

where  I  denotes  the  total  number  of  pulses  transmitted  per 
burst,  Tr  is  the  pulse-repetition  interval  (PRI),  C  denotes  the 
transmitted  code  matrix  that  is  being  optimized  and  s(t)  is  a 
fixed  unit  energy  “template”  pulse  emitted  by  each  transmitter. 

A.  Data  Model 

Over  Nb  range  bins,  the  received  signal  can  be  represented 
as  a  NrI  x  Nb  matrix  Y  where  each  column  is  a  length- AVr  / 
vector  corresponding  to  a  single  range  bin  and  is  written  as 

y  =  t  +  q  +  w,  (2) 
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where  {t,  q}  denote  the  signal  contributed  by  the  target  and 
clutter,  and  w  represents  additive  Gaussian  noise  with  known 
covariance  matrix  R„, . 

Traditional  models  for  the  target  generally  assume  a  point 
target  with  freespace  line-of-sight  (LOS)  propagation;  simi¬ 
larly,  clutter  is  modeled  as  a  superposition  of  freespace  LOS 
propagation  from  an  assortment  of  clutter  sources,  e.g.,  [12]. 
For  simplicity,  we  consider  a  two-dimensional  radar  configu¬ 
ration  in  the  xy  plane,  with  the  radar  signal  ray  path  forming 
an  angle  9  with  respect  to  the  x-axis. 

In  contrast  to  the  freespace  case,  propagation  in  inhomo¬ 
geneous  media  suffers  from  the  effect  of  refractive  index 
fluctuations,  leading  to  scintillations  in  the  signal  amplitude 
and  phase.  We  focus  on  the  phase  scintillation  in  this  paper 
as  this  effect  is  the  dominant  influence  on  the  temporal 
coherence  (Doppler  spread)  and  spatial  coherence  (angle-of- 
arrival  spread)  of  the  radar  signal.  The  phase  of  a  signal  at 
a  given  point  in  time  is  equal  to  the  integration  of  its  spatial 
rate  of  change  over  the  ray  path: 

</>  =  f  k(r)ds,  (3) 

j  So 

where  k{ r)  is  the  radar  signal  wavenumber,  ds  is  an  element 
of  arc  length,  and  (so,si)  are  the  ray  path  end  points.  If  we 
suppose  there  are  no  zero-order  plasma  density  gradients  in 
our  assumed  two-dimensional  plasma  volume,  then  the  rays 
form  straight  lines  that  can  be  taken  as  radial  with  respect 
to  the  radar  in  a  system  of  polar  coordinates  (p.  9).  Let  us 
consider  a  first-order  Taylor  series  pertubation  of  k( r)  with 
respect  to  plasma  density  n,  such  that  we  can  write  the  phase 
perturbation  as 


where  n\  is  the  plasma  density  perturbation,  and  (po,  p\)  is  the 
interval  of  integration  in  the  radial  coordinate.  The  dispersion 
relation  for  an  unmagnetized  plasma  is  given  by 

2  2  7  2  ,  rc\ 

UJ  =  C  k  - 1 - ,  (5) 

e0me 

where  ui  is  the  carrier  frequency,  c  is  the  speed  of  light,  qe  is 
the  charge  on  an  electron,  eo  is  the  permittivity  of  free  space, 
and  me  is  the  mass  of  an  electron.  Using  the  dispersion  relation 
to  compute  dk{r)/dn,  we  arrive  at 

rPi 

</>i=~reX  ni(r)  dp,  (6) 

J  PO 

where  re  =  q2/ (AireQmec2)  =  2.8  x  10”15  m  is  the  classical 
electron  radius,  and  A  is  the  radar  wavelength  in  the  plasma. 
The  spatial-temporal  autocorrelation  of  <j> i  in  the  assumed  two- 
dimensional  geometry  is  given  by 

R<t>±  (X,Y,T)  = 

(reX)2  f  f  Rni(X  +  x  —  x',Y  +  y  —  y',T)  dp  dp'.  (7) 

d  Pa  J  Pa 


After  Fourier  transforms  we  have 


(^77  1  Ky  ,  12)  -  (reA)  Sni  (nx:  Ky  ,  fl) 

x  rpi  rPi  ei(Kx  cosfl+Kj,  sinff)(p-p')  ^  ^  / 
dpo  dp0  "  " 

=  2nL(reX)2S(nx  cos  9  +  Ky  sin 9)Sni (kx,kv,  fl),  (8) 


where  L  =  pi  —  po  is  the  path  length,  and  it  is  observed  that 
the  Dirac  delta  function  assures  that  ny  =  —Kxcot  9.  For  a 
two-dimensional  turbulent  plasma,  the  Kolmogorov  turbulence 
spectrum  has  a  k~8^3  dependence,  which  we  approximate  as 


.-3. 


^n1  (ft’X  j  tty ;  — 


Ait2ko  (n\ )  5(Q) 


("g- 


2)3/2  ’ 


(9) 


where  kq  is  the  inverse  scale  length  of  the  largest  turbulent 
eddies  and  we  have  assumed  no  temporal  structure.  The 
spectrum  is  normalized  as  follows: 


<»?>  = 


(2tt)3 


Sni  ( kx ,  Ky ,  f l)  dnx  dny  dfl. 


(10) 


If  the  plasma  is  drifting,  then  temporal  structure  is  created  by 
the  movement  of  spatial  structure  past  a  point  in  space,  namely 


Air2 Kq  (nl)  6(Q  -  Kxvdx  -  Kyvdy ) 

(Kg  +  «£  +  «3)3/2  ■ 


(ii) 


By  inserting  this  expression  for  Sni(nx,  Ky,fl)  into  (8)  and 
recalling  that  Ky  =  —Kxcot  9,  we  can  integrate  out  the  kv 
coordinate  so  that  the  phase  spectrum  becomes 


S(p i  (kx  ,  D)  — 

8ir3KoLcsc9(reX)2  (n3)  6  (D  —  Kx(vdx  —  vdy  cot  9)) 
(Kq  +  k2  esc2  9)3/2 


(12) 


Given  the  phase  spectrum  in  (12),  the  data  cube  is  then 
crated  by  forming  a  two-dimensional  sensor-pulse  phase  signal 
at  each  range  bin  by  filtering  white  noise  through  H  = 
IFFT(yjS^).  The  resulting  data  cube  is  then  convolved  with 
the  ambiguity  function  of  the  transmitted  signal. 


B.  Target  Contribution 

The  target  is  assumed  to  lie  in  a  single  far-held  range  bin 
and  that  it  can  be  modeled  as  a  point  source  with  amplitude 
at.  Further,  the  target  is  located  at  an  azimuth  angle  9t  with 
respect  to  the  receive  array  and  Doppler  frequency  ilt.  In  this 
application,  the  key  difference  from  previous  approaches  is 
that  the  propagation  through  plasma  media  generates  random- 
phase  perturbations  in  the  signal  contributed  by  the  target  at 
the  receiver.  These  perturbations  must  be  accounted  for  in  the 
design  of  the  transmit  waveform  and  adaptive  receiver. 

Let  <p\n  represent  the  random  phase  perturbation  at  the  n- th 
(of  Nr)  receiver  and  the  v'-th  (of  J)  pulse.  Based  on  the  target 
model  in  Section  II-A,  after  matched  filtering  to  the  transmit 
pulse,  the  (m)-th  entry  of  t  can  be  written  as: 
nt 

tin  =  at  X]  Cmiej^ej2M^ej2^m+n\  (13) 

m=  1 

where  7  =  dx/dn,  fs  =  9r  cos(0t)/A  is  the  spatial  frequency, 
and  Tr  is  the  Pulse  Repetition  Interval  (PRI). 
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For  convenience,  we  define  the  following  matrices  and 
vectors: 

•  the  length -INt  vector  c  =  VEC(C)  and  ( NrINt )  x 
(NrINt)  augmented  transmit  code  matrix  as  C  = 

INr  <g )  C; 

•  the  Doppler  vector  a<>  as  the  length-/  vector  [ao],  = 
exp(j27riOtTr).  Further,  define  the  I  x  I  diagonal  matrix 
An  =  diag(an); 

•  based  on  these  definitions,  define  two  forms  of  a  spatial- 
temporal  matrix 

0  =  a^(g)An®aT,  (14) 

©  =  diag(aR)  ®  An  ®  (aT)',  (15) 

with  a t  and  a r  being  the  transmit  and  receive  steering 
vectors,  respectively. 

•  finally,  define  length-/ A/j  vector  pt  such  that  \i>t\i+i,r  = 
<t\n,  the  vectorized  form  of  the  phase  perturbations  and 
the  INR  x  IN n  diagonal  phase  perturbation  matrix  Pt 
as  [P t\in  =  diag(pt). 

Based  on  these  definitions,  the  target-generated  signal  in  (13) 
can  be  re-written  in  the  two  following  forms: 

t  =  at0Cp  t,  (16) 

t  =  aiPt0c.  (17) 

C.  Clutter  Contribution 

As  seen  in  Section  II-A,  the  clutter  model  is  similar  to 
the  target  model  with  multiple  clutter  rays  contributing  to  the 
overall  signal.  Using  (12),  we  know  that  the  clutter-generated 
phase  perturbations  are  both  spatially  and  temporally  corre¬ 
lated.  Inspired  by  the  signal  model  proposed  for  the  target, 
the  clutter-contributed  signal  can  be  also  written  as 
Ne 

q  =  ^09(0„)Cp9(0w),  (18) 

V  —  l 

where  Ng  denotes  the  number  of  azimuth  paths  over  which 
the  clutter  signal  is  received,  and  Qq(9v)  =  diag(aft(0„))  8 
An,g(0„) <g>  (aT(9v))'  and  pq(0v)  defined  similarly  to  the  cor¬ 
responding  terms  for  the  target  signal  but  for  clutter  azimuth 
angle  0V.  Similarly,  we  have  A n,q(9v)  =  exp(j2niSl(9v)Tr) 
where  f2(0„)  denotes  the  mean  Doppler  shift  of  the  clutter 
patch  at  azimuth  angle  9V. 

We  assume  that  different  azimuth  angles  contribute  phase 
perturbations  that  are  statistically  independent.  Specifically, 
given  4)1  (0V)  and  (pq,.-, (9V)  as  the  perturbations  generated 
by  the  u-th  and  u'-th  path,  respectively,  E{(f>L(9v){(j)q)  p 
(9V>)}  =  0,  Vu  ^  v' .  Based  on  this  assumption,  we  focus  for 
now  on  a  single  ray  path  at  azimuth  angle  9.  Define,  for  an 
array  along  the  x-axis,  the  spatial-temporal  correlation  of  the 
phase  perturbation  as 

R(Ax,  At)  =  E  {4>{x,  t)(f>*(x  +  Ax ,  t  +  A/)}  ,  (19) 

where  Ax  and  At  are  the  spatial  and  temporal  lag,  respec¬ 
tively,  and  E{-}  denotes  expectation.  The  phase  autocorrela¬ 
tion  matrix  RA  is  a  (NppI)x(NpiI)  matrix  whose  ( i+nl ),  (i'+ 
n'l)- th  entry  is  as  R Xi+ni)(i'+n'i)  =  R(\n-n’\dR,  \i-i'\Tr). 


Note  that  the  clutter  contribution  in  (18)  is  a  nonlinear 
function  of  the  phase  and,  therefore,  further  analysis  is  needed 
to  obtain  the  autocorrelation  matrix  of  received  signals.  We 
use  the  fact  that  the  phase  perturbations  &)n  is  as  a  two- 
dimensional  Gaussian  random  process. 

Lemma  1.  Let  zin  =  [pq]i+np  =  exp(j<^)fn)  represent 
the  clutter  phase  term  in  (18).  The  clutter  phase  auto¬ 
correlation  matrix  is  a  (INpp)  x  ( /JVr )  matrix  R1^  = 

-  ^{P?)}(P<?  ~  MpJ)1'}  with 

T>  4> 

E  =  E  ZinZ* i'n' 

=  exp  (~a2q)  (exp  -  l)  ,  (20) 

where  z  =  E{z}),  a 2  is  the  clutter  variance,  and  = 

Rp(|n  -  n'\dR,  \i  -  i’\Tr). 

The  autocorrelation  matrix  of  the  clutter  signal  given  in  (18) 
for  a  clutter  ray  from  azimuth  angle  9  can  be  then  written  as 

R,(0)  =  $(0)R0($(0))t,  (21) 

with  $(0)  =  09(0) C. 

Finally,  for  the  more  practical  case  that  the  clutter  is  re¬ 
ceived  over  Ng  >  1  rays  and  using  the  statistical  independence 
of  different  rays,  we  have 

Ne 

R9  =  ^R?(0i;),  (22) 

V—l 

where  each  individual  Rg(0„)  is  given  by  (21). 

III.  Joint  Optimization  of 
Signal-to-Interference-plus-Noise  Ratio 

At  the  receive  array  the  signals  received  within  each  range 
bin  are  adaptively  combined  to  maximize  the  SINR.  Detection 
or  localization  decisions  are  then  made  on  the  filter  output. 
Denoting  as  h  the  length- NRI  vector  of  adaptive  weights,  the 
SINR  is  given  by 

E  { |  h1 1 1 2  } 

SINR  =  „  f„  t  — '  *  ,  ,91  .  (23) 

/?{|htq|2}  +  i?{|hfw|2} 

Note  that  the  target  and,  importantly,  clutter  signals,  t  and  q, 
are  functions  of  the  transmit  code  vector  c.  We  are  interested 
in  jointly  designing  the  unknown  adaptive  filter  weight  vector 
h  and  transmit  signal  c  under  a  constraint  on  the  total  trans¬ 
mitted  power.  This  optimization  problem  can  be  formulated  as 
follows: 

{c°,h°}  arg maxh.c  SINR, 
such  that  llcl|2  5s  1- 

The  cost  function  in  (24)  is  nonlinear  and  non-convex  with 
respect  to  the  joint  optimization  variables  [cl  hi] I.  Globally 
optimal  solutions  are,  therefore,  impossible  to  come  by  and 
here  we  propose  an  iterative  approach  in  which  the  linear 
weight  vectors  h  and  input  signal  vector  c  are  obtained 
individually. 
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A.  Optimization  of  Adaptive  Weights 

In  optimizing  for  the  adaptive  weights,  we  assume  the 
transmit  waveform,  defined  by  c,  is  fixed.  For  convenience,  the 
formulation  in  (16)  is  used  to  optimize  the  adaptive  weights 
h.  The  SINR  can  be  then  written  as  follows: 


SINR  = 


H2htpc,th 

htPc,9,™h 


(25) 


where  the  subscripts  t  and  q  refer  to  the  target  and  interfer¬ 
ence  (clutter  and  noise)  respectively,  and  pCyt  and  pc,q,w  are 
two  known  matrices  that  are  functions  of  target  and  clutter 
parameters.  Here, 


Pet  =  0tCR?(C)t(0t)t 


(26) 


Pc,q,w  Rc  T  Ru 

No 


=  e,(0„)CR^„)(C)t(e,)t(0„)  +  RW1  (27) 


V  =  1 


where  Qq(9v)  and  R^(9V)  are  calculated  using  (15)  and  (20), 
respectively,  for  each  clutter  azimuth  angle  6V.  In  addition, 
R„,  denotes  the  covariance  matrix  of  the  additive  noise. 

The  optimal  weight  vector  (h°),  for  a  given  transmit  code 
c,  can  be  then  found  by  maximizing  the  cost  function  in  (25). 

A  more  stable  approach  is  to  use  the  Cholesky  factorization 
of  p  .  Defining  Lc,q,w  as  the  Cholesky  factor  of  pc 
where  pcqw  =  Lc,q,w(L  c,q,w)>  the  following  lemma  presents 
a  framework  for  finding  the  optimal  solution  of  the  liner 
weights: 


Lemma  2.  The  optimal  solution  of  the  problem  in  (25)  is: 


where  'P(A)  corresponds  to  the  principal  component  of  matrix 

A. 


B.  Waveform  Optimization 

The  second  step  within  the  iterative  algorithm  reverses  this 
to  find  the  optimal  waveform  given  the  adaptive  weights  h.  In 
this  case,  using  the  alternate  signal  model  in  (17),  the  SINR 
can  be  written  as 


SINR  =  |  at  | 


c^rytc 

ctfl,ljCc  +  (h)tRUJh 


(29) 


where  f 2^*  and  0/,  c  are  known-structure  matrices  function 
of  target  and  clutter  parameters,  respectively.  The  waveform 
design  can  be  now  reformulated  as  the  following  optimization 
problem: 

{c°}  =  arg  maxc  SINR, 
such  that  =  1 1 c| | 2  <  1. 


It  was  shown  in  [11]  that,  within  a  scale  factor,  the  above  op¬ 
timization  problem  leads  to  the  same  solution  as  the  following 
unconstrained  SINR-maximization 


{c°}  =  arg  max 

C 


_ c*Slh,tC _ 

ctf2/,jCc  +  (h)fRtoh.cfc 


(31) 


If  c°  is  a  solution  obtained  for  this  unconstrained  problem,  the 
final  solution  is  c°  =  c°/||c°||. 

The  cost  function  in  (31)  can  be  simplified  into  the  follow¬ 
ing  form: 


SINR  =  |at|2 


C^h,tC 

^h,c,wC 


(32) 


with  flh,c,w  =  nh.c+ttytR^hljv*/  and  solution  obtained  in 
a  manner  similar  to  the  optimization  problem  for  the  known- 
waveform  case.  Define  V h,c,w  as  the  Cholesky  factor  of 
f lh,c,w  Now,  using  lemma  2  the  optimal  waveform  is  found 
as  follows: 


c,  = 


V”1 

h,c,w 


V 


[V  h,c,w] 


h,t 


v 


-i 


it 


(33) 


The  final  waveform  is  then  obtained  by  normalizing  the 
solution  obtained. 


C.  Overall  Joint  Optimization  Algorithm 

The  overall  optimization  algorithm  is  as  follows: 

1 )  Given  the  phase  PSD  for  the  clutter  and  target,  calculate 
the  autocorrelation  matrices  R f  and  R %(6V)  where  v  = 
1,  •  •  •  ,N$  and  Ng  denotes  the  number  of  clutter  rays. 

2)  Based  on  the  target  and  clutter  parameters,  compute  the 
spatial-temporal  matrices  0t,  <dt,  Qq(9v)  and  Qq(9v). 

3)  Initialize  matrices  pct  =  0tCRf (C)^(0t)^  and 

Pc,q,w  =  Evli  6,(0„)’cR|(fl„)(C)t(6,)t(eu)  +  rw. 

4)  Find  the  Cholesky  factor  of  pcqw  as  hc,q,w 

5)  Calculate  the  optimal  filter  weight  using  (33)  by  h°  = 


V 


[Lc^uj]  Pc  t  [-h,Clg,u>] 

6)  Form  the  new  matrix  H°  =  diag(h°). 

7)  Initialize  new  matrices  flh,t  =  0^HoRf (H°)^0t 
and  SlhiC'W  =  Ef=i0j(^)H°R^(^)(H°)t09(^)  + 
(h°)tRu,h°IjvJj/. 

8)  Let  V h,c,w  represent  the  Cholesky  factor  of  flh,c,w- 

9)  Using  (33),  set 


cr  = 


vr1 

ri,c,w 


it 


v 


Vr1 

h,c,w 


as  the 


[Vh,c,iu]  flh,t 

optimal  solution  for  the  modified  optimization  problem. 

10)  Calculate  c°  =  c°/||c°||  as  the  optimal  waveform. 

11)  Repeat  steps  2)-10)  until  convergence,  defined  as  a 
negligible  change  in  the  resulting  SINR. 


IV.  Numerical  Results 

We  consider  a  scenario  with  a  single  target  at  an  azimuth 
of  70°  and  Doppler  of  1/6  of  the  pulse  repetition  frequency. 
There  are  4  transmit  antennas  and  4  receive  antennas,  placed  at 
half-wavelength  spacing.  The  waveform  consists  of  a  coherent 
burst  of  I  =  36  pulses.  In  simulating  and  modeling  distributed 
clutter,  an  azimuth  angle  step  of  59q  =  1°  is  used.  We  use 
the  clutter-to-noise  Ratio  (CNR)  to  characterize  the  power  of 
the  incoming  clutter  in  the  receiver.  The  additive  noise  is  a 
circularly  symmetric  complex  Gaussian  random  variable  with 
unit  variance. 

The  performance  of  the  waveform  design  framework,  and 
the  benefits  of  waveform  design,  is  now  evaluated  by  analysing 
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et  =  70,  CNR  =  50  (dB) 


0(  =  70,  a  =  20  (dB) 


Fig.  2.  ROC  curves  for  optimal  waveform  and  Hadamard  code  designs. 


Fig.  1.  The  beampattern  of  the  optimal  waveform  for  the  single-ray  clutter 
and  two  different  distributed  clutter  cases. 


the  beampattern  of  the  optimal  code.  Figure  1  shows  the 
performance  of  the  joint  waveform-receiver  adaptive  process 
in  the  practical  case  of  distributed  clutter.  The  figure  plots 
the  adapted  beampattern  for  a  single  clutter  ray  ( 8q  =  30°) 
and  two  cases  of  distributed  clutter  (30°  <  9q  <  40°  and 
30°  <  9q  <  50°).  As  before,  the  target  is  at  9t  =  70°.  Clearly, 
the  algorithm  is  able  to  adapt  to  the  clutter  by  placing  broad 
null  in  the  required  angular  region.  The  penalty  for  placing  a 
broad  null  is  also  evident  in  the  reduced  gain  on  the  target. 

The  ability  to  place  a  null  in  direction  of  the  clutter 
improves  detection  performance.  We  now  compare  the  detec¬ 
tion  performances  of  the  waveform  optimization  and  fixed- 
waveform  schemes.  In  each  case,  assuming  Gaussian  signals, 
the  probability  of  detection  (Pd)  is  related  to  probability  of 
false  alarm  (Pfa)  as 


Pd  =  1  —  F, 


*2 


^  1  +  SINR  )  ’ 


(34) 


where  Fx2  (•)  denotes  the  %2  Cumulative  Distribution  Function 
(CDF)  with  two  degrees  of  freedom. 

Figure  2  plots  the  radar  receiver  operative  characteristic 
(ROC)  curve,  comparing  the  performance  of  the  optimized 
waveform  with  that  of  using  a  Hadamard  code.  The  clutter 
is  distributed  over  the  span  of  30°  <  9q  <  50°.  Again, 
both  schemes  use  adaptivity  on  receive  and  the  improvements 
are  due  to  optimizing  the  transmit  waveform.  The  huge  im¬ 
provements  in  detection  performance  are  evident,  especially 
at  the  lower  probabilities  of  false  alarm.  For  example,  for 
Pfa  =  10”6,  the  optimal  waveform  provides  a  detection  rate 
of  96.5%  while  the  Hadamard  code  only  provides  a  58.5% 
detection  rate. 

Figure  3  plots  the  detection  probabilities  as  a  function  of 
target  amplitude  (a),  again  comparing  the  performance  of  the 
optimal  waveform  and  Hadamard  code.  For  comparison,  the 


Fig.  3.  Probability  of  detection  versus  target  amplitude  for  optimal  waveform 
and  Hadamard  code  designs. 


detection  curves  were  obtained  for  two  cases;  one  with  single¬ 
ray  clutter  assumption  and  the  other  with  distributed  clutter 
where  30°  <9q<  50°.  In  both  cases,  the  CNR  is  set  to  50dB. 
The  improved  detection  performance  is  again  clear,  especially 
in  the  case  of  distributed  clutter.  The  figure  also  illustrates  the 
performance  penalty  in  nulling  distributed  clutter.  The  results 
in  Fig.  3  show  that  using  a  fixed  transmission  code  is  more 
sensitive  to  the  structure  of  the  clutter  whereas  optimizing  the 
waveform  makes  the  detection  performance  more  robust.  As 
expected,  at  very  high,  likely  impractically  high,  SNR  levels 
the  performance  gap  is  eliminated. 

Our  final  set  of  results  focus  on  how  the  optimal  waveform 
design  improves  Doppler  frequency  estimates  as  compared 
to  non-optimal  waveforms.  In  this  example,  the  transmitters 
send  /  =  36  pulses  per  coherent  pulse  interval  leading  to 
36  Doppler  resolution  cells.  Doppler  data  are  now  generated 
for  the  optimal  waveform,  Hadamard  code  and  random  code. 
Assuming  CNR  =  30dB,  a  =  14dB  and  the  distributed  clutter 
case  with  30°  <  9q  <  50°,  Figure  4  plots  the  received 
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power  in  different  Doppler  cells  for  all  three  transmission 
schemes.  The  improved  Doppler  discrimination  using  the 
optimal  waveform  is  clear.  The  Doppler  power  spectrum  for 
the  optimal  waveform  has  a  clear  peak  at  the  target  Doppler 
frequency  (flt  =  PRF/6).  Although  the  Hadamard  code  also 
shows  maximum  power  around  the  same  Doppler  cell  that  was 
caught  by  the  optimal  waveform,  its  power  is  still  4dB  lower 
than  the  maximum  power  achieved  by  the  optimal  code.  The 
Hadamard  code  also  leads  to  spurious  Doppler  peaks.  Finally, 
unlike  two  aforementioned  schemes,  the  random  code  fails 
to  detect  the  true  Doppler  frequency  as  the  maximum  power 
occurs  at  incorrect  Doppler  frequencies. 
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Fig.  4.  The  received  power  in  different  Doppler  cells  for  the  optimal 
waveform,  Hadamard  code  and  random  signal. 


V.  Conclusions 

This  paper  has  proposed  a  technique  for  the  joint  design  of 
the  waveforms  and  adaptive  processor  in  a  radar  system.  Using 
the  model  for  the  phase  perturbations  to  derive  related  signal 
and  clutter  models,  the  transmitted  waveform  and  matched- 
filter  weights  are  optimized  to  maximize  the  output  SINR  of 
the  system.  Importantly,  and  unlike  some  previous  approaches, 
the  impact  of  the  transmitted  waveform  is  accounted  for  in 
modeling  the  clutter.  Finally,  we  present  simulation  results 
to  illustrate  the  performance  gains  of  the  optimal  waveform 
design  over  other  non-optimal  waveform  techniques. 
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